Metal— Insulator Transition in 2D: Experimental Test of the Two-Parameter Scaling 
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We report a detailed scaling analysis of resistivity p(T, n) measured for several high-mobility 
2D electron systems in the vicinity of the 2D metal-insulator transition. We analyzed the data 
using the two parameter scaling approach and general scaling ideas. This enables us to determine 
the critical electron density, two critical indices, and temperature dependence for the separatrix in 
the self-consistent manner. In addition, we reconstruct the empirical scaling function describing a 
two-parameter surface which fits well the p(T, n) data. 

PACS numbers: 71.30.+h, 73.43.Nq, 71.27.+a 



The experimental discovery [l|, 0, 0] of the metal- 
insulator transition (MIT) in the two-dimensional (2D) 
electron system of high mobility silicon metal-oxide- 
semiconductor (Si-MOS) field-effect transistors still calls 
for a theoretical explanation. Perhaps, the most promis- 
ing framework is provided by the microscopic theory, ini- 
tially developed by Finkelstein, which combines disorder 
and strong electron-electron (e-e) interaction [j]. Re- 
cently, Punnoose and Finkelstein [5(] have shown the pos- 
sibility of the MIT for a special model of the 2D electron 
system with the infinite number of valleys. The key fea- 
ture of the theory is the existence of the two-parameter 
scaling near the criticality. One scaling variable is gov- 
erned by disorder, and the other is determined by e-e 
interaction. In accord with the theoretical prediction, 
the renormalization of both resistivity and e-e interaction 
with temperature has been experimentally demonstrated 

The main objective of the present Letter is to test 
experimentally the two-parameter scaling. We present 
our results for the detailed analysis of the temperature 
(T) and electron density (n) dependence of resistivity 
(p) measured on five Si-MOS samples. From this anal- 
ysis, we identify the two scaling variables together with 
the corresponding exponents (see Table [TJ) and determine 
for the first time the scaling function for resistivity in 
the metallic region in the vicinity of the MIT. We have 
found that the two-parameter scaling provides a natural 
explanation for the following set of experimental obser- 
vations in the vicinity of the MIT: (i) the separatrix, 
p c (T) = p(T,n c ), has a monotonic power-law tempera- 
ture dependence, where n c denotes the critical density; 
(ii) a generic p(T, n) curve on the metallic side (n > n c ) 
is nonmonotonic in T, and has a maximum and an in- 
flection point; (iii) for n — n c , the maximum and inflec- 
tion points merge at zero temperature; (iv) the normal- 
ized p(T,n)/ p c (T) data in the vicinity of the MIT obey 
the mirror reflection symmetry that is fulfilled in much 
wider ranges of n — n c and T, as compared to the sym- 
metry previously reported under the assumption of the 



T-independent separatrix [8j. 

For the measurements we have selected five repre- 
sentative Si-MOS samples of the rectangular geometry 
from four different (OOl)-Si wafers. Their peak mobil- 
ities are listed in Table HI We have used four-terminal 
ac-technique at the 5 Hz frequency. The source-drain 
current was chosen low enough, e.g., I = 10 nA for 
T = 1.3 K, in order to avoid electron overheating. Most 
of studies were performed in the range 1.3 — 4.2 K, be- 
cause for these temperatures (i) the critical region of den- 
sities | An | = | ri — n c \ becomes sufficiently wide as will be 
shown below, and (ii) p(T) near the criticality exhibits 
a well-pronounced maximum (Fig. [1]) , which provides an 
additional tool for comparison with the theory. For two 
samples, Sil5 and Si62, we have performed the measure- 
ments down to 0.3 K; for sample Si43, we have extended 
the temperature range up to 1 — 38 K in order to illu- 
minate the overall trend of the resistivity maximum and 
the scaling behavior of p(T, n). 

Figure [T] shows a typical o(T)-dependence for two sam- 
ples, Si2 (Si43), in the critical regime, in the density range 
from 0.672 to 1.12 (0.707 to 1.49), in units of 10 11 cm" 2 
[9j. As n decreases, the p(T) behavior crosses over from 
metallic (n > n c ) to insulating(n < n c ) one. Initially 
and during about a decade 0, Hf, HI , the data were an- 
alyzed within the framework of the one-parameter scal- 
ing (OPS) theory, in which n is a single driving param- 
eter. Within this approach the separatrix p c (T) sepa- 
rating the metallic from insulating domains should be 
T-independent, i.e, p c (T) = po. As a result, in the 
vicinity of n c , the p(T,ri) data should obey the reflec- 
tion symmetry [llj: p(T,n c — An) = p^/p(T,n c + An) 
when An <C n c . This obviously disagrees with the data 
shown in Fig. [TJ i.e., the OPS theory cannot adequately 
describe the data within the whole T-range. In view of 
this well-known problem 12j , it is a common practice to 
use the OPS approach for the analysis of the data in a 
truncated temperature range. 

In theory physics at T = is described by the 
two coupled renormalization group (RG) equations for 
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conductance a (in units e 2 /h) and interaction amplitude 
72 = -Fg/(1 + F$) where F CT denotes the standard 
Fermi-liquid parameter in the triplet channel 0, EH : 

^ =A 7 (o-,7 2 ), ^ = /? 72 (< J ,72)- (1) 

Here 77 = InL/Z, L stands for the running RG length- 
scale (sample size), and I is the microscopic lengthscale 
(elastic mean free path), at which the RG starts. In 
the microscopic theory [J], the temperature appears as 
Tz which has a dimensionality of L~ 2 , and the param- 
eter z varies with L according to dhiz/drj = 72(17,72). 
The RG equations should be supplemented with the ini- 
tial conditions at r\ — 0: cr(0) = ct, 72(0) = 72, and 
z(0) = z, where (7,72 and z are sample dependent. Fol- 
lowing the scaling ideas we expect the measured 
resistivity p — p(Tz, a, 72, 1) to be independent of micro- 
scopic details: dp/dl = 0. Applying the standard analysis 
of the Callan-Symanzyk equation [14| . we find that the 
measured resistance is a function of two scaling variables 
X and Y, p = 1Z(X, Y), where TZ is a regular function. 
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FIG. 1: Temperature dependences p(T): (a) for sample Si2; 
the densities increase in steps of 0.224 starting from 6.72 (top- 
most curve), (b) for sample Si43; the densities are (from top 
to bottom) 7.16, 7.34, 7.52, 7.70, 7.88, 8.06, 8.51, 8.96, 9.41, 
9.86, 10.31, 10.76, 11.6, 12.5, 13.4, 14.9. The inset shows T max 
vs electron density: squares are the data, and line is the the- 
oretical curve (see below). Density values are quoted in units 
ofl0 10 cm- 2 . 

The MIT implies the existence of a fixed point (FP) 
determined by the conditions: /3 CT (cr c , 72) = /3 72 (<r c , 7f) = 



0. Linearizing Eqs. |T]) near the fixed point, one finds 

X = {T/T )- K (n - n c )/n c , Y = (T/T t )<. (2) 

Here k — p/(2v) and Q = —py/2, where the correlation 
length exponent v and the irrelevant exponent y < are 
given by the eigenvalues of the linearized Eqs. ([T]) 
The exponent p = 2/[2 + 7z(ct c ,7|)] governs the T- 
behavior of the specific heat at the fixed point; Tq^i are 
sample-dependent energy scales. Temperature To corre- 
sponds to the bandwidth of the available states partici- 
pating in the transport and is of the order of the elastic 
scattering rate fi,/2fcsr. Temperature T\ determines the 
quantum critical region at n — n c [i~5| . Deriving Eq. l[2]). 
we expressed a and 72 in terms of the experimentally 
measurable parameters p and n. 
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FIG. 2: The density dependences of In p(n) for vari- 
ous temperatures (from the top to the bottom, in K): 
(a) 1, 1.23, 1.43, 1.6, 1.8, 2.1, 2.35, 2.64, 2.95, 3.31, 3.71, 4, 4.33 
for sample Si43, (b) 1.4, 1.6, 1.8, 2.2, 2.5, 3, 3.6, 4.2 for sample 
Si2. The density is given in units of 10 11 cm -2 . The data 
are shown in a limited density range, X < 0.7. The insets 
illustrate the procedure of extracting k and To using Eq. I|3). 
The bars display the range of intersections for all curves. 

At present, the explicit form of 1Z(X, Y) is unknown. 
In general, even in the vicinity of the FP, X and Y are 
not necessarily small. However, from the common-sense 
arguments we require that at T — > the system (i) is 
an insulator for n < n c , i.e., TZ(— 00, 0) = 00; (ii) is a 
metal for n > n c , i.e., 7£(+oc,0) = 0; and (iii) is at the 
quantum critical point for n = n c , i.e., 7Z(0,Q) = p® c . 
Following arguments of Ref . [llj , we assume that 

p(T,n) = K(X,Y)=p° c e- x (l~Y), \X\, Y < 1. (3) 
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This is consistent with all the requirements listed above. 
In particular, such a form of 1Z(X, Y) provides the reflec- 
tion symmetry for the normalized resistance p(T)/p c (T). 

In the vicinity of n c , on the metallic side, n > n c , 
the p(T, n) data exhibit a maximum at a certain density- 
dependent temperature T max (see Fig. [1]) . The inset to 
Fig. [T] shows that T max vanishes as n — > n c , in agreement 
with Eq. ([3]), which predicts a monotonic p c {T) depen- 
dence. On the insulating side, for n lower but close to 
ri c , the p(T,n) curves have a positive curvature. There- 
fore, the neighboring curves on the metallic side should 
also have a positive curvature at sufficiently high tem- 
peratures. This expectation agrees with the experiment 
(Fig. [TJ. It follows from the above results that each of 
the nonmonotonic p(T, n) curves must have an inflection 
point. Its existence restricts the exponent £ to be less or 
equal to unity, C ^ 1- As a result, not only T max but also 
the inflection point temperature Tj n f vanish at n = n c . 

In order to test the theoretical prediction, we have 
adopted a procedure employed for studying the scal- 
ing behavior near the criticality in the quantum Hall 
regime [lj| . We plot in Fig. [5] the In p(T, n) data for 
various temperatures versus n for X < 1. Two features 
are clearly seen in Fig. [2] Firstly, as T decreases, the 
region of intersection for the lnp(T, n) curves shrinks. 
This enables one to determine an approximate n c value 
as an intersection point in the T = limit. Secondly, 
in the vicinity of n c (for small X), the lnp(T, n) curves 
are linear in n in accordance with Eq. ([3]). In the insets 
to Fig. [2j we plot the logarithmic slopes of the In p(T, n) 
curves versus InT. The resulting curves appear to be 
linear; this allows us to extract both k and To values for 
the five studied samples. The results are summarized in 
Table H 

Equation ([3]) predicts that, for X <C 1, the p(T,n) 
data normalized by exp(— X) should collapse onto a sin- 
gle curve. In order to test this theoretical prediction, 
we plot the normalized p*{T) = p(T,n) exp(X) data in 
Fig. [3J We stress that no adjustable parameters are used 
in this scaling procedure. The critical density n c , and the 
To and k values were determined at the previous stage. 
As is shown in Fig. [3l the normalized p*(T) data do scale 
with high accuracy, 0.8%, on the p(T,n c ) curve in the 
range X < 0.5. The deterioration of the scaling quality 
for greater X values limits the range of temperatures and 
densities used for the scaling procedure (Fig. [3]). Pro- 
vided that Eq. (J3j) adequately describes the separatrix 
p(T,n c ), the scaling procedure enables us to determine 
the critical n c value more precisely. For this purpose, we 
require the low temperature region (Y <C 1) of the scaled 
curves to have the curvature of the constant sign over 
the whole temperature range. Indeed, we found that, 
by adjusting the preliminary determined n c value within 
only 2%, the inflection point in the p(T, n c )-data could 
be easily eliminated for all five samples. 

The bunch of the scaled curves in Fig. [3] is fitted with 



the theoretical expression ([3]) for the separatrix p(T, n c ) , 
by using the remaining scaling parameters £, T\ and p® c 
as adjustable parameters. The result of fitting for sam- 
ple Si2 is shown in the inset to Fig. [3] The values of the 
above parameters extracted from the best fit are given in 
TableUJfor all of the five samples. The validity of Eq. © 
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FIG. 3: Scaling of the p*(T)/p° data for densities n = 
8.74, 8.96, 9.18, 9.41, 9.63 x 10 10 cm" 2 . The lower inset demon- 
strates the quality of the scaling for large X ~ 0.5. The upper 
inset shows a fit of the separatrix using Eq. (|3} with three pa- 
rameters £, Ti and p°. Sample Si2 

is further supported by the nonmonotonic temperature 
behavior of p(T, n) (Fig. [JJ). Using the extracted scaling 
parameters (Table[TJ), we can readily find T max (n) depen- 
dences from Eq. The results are in a good agreement 
with the experimental data for all of the five samples. A 
typical example of comparison is shown in the inset to 
Fig. [TJ We note that the data are processed only for 
T < (4 - 5) K to ensure the applicability of Eq. (J3J. 

TABLE I: The relevant parameters for the studied samples. 
Ppeak [m 2 /Vs] is the peak mobility at T = 0.3 K, p° c is in 
units of h/e 2 , To and Ti are in K, and n c is in 10 11 cm~ 2 . 
The inverse transport scattering time fi/2fcsr [K] is estimated 
from p c (T = 4.2 K). The errors for To, k, n c are about 2%; 
the errors for T\, C, and p° are given in the last column. 



sample 


/^pcak 


h/2k B T 


n c 


To 


Ti 


K 


c 


P° c 


% 


Sil5 


4.1 


13 


0.86 


12 


15 


0.82 


0.35 


4.17 


30 


Si62 


3.6 


26 


0.94 


22 


11 


0.74 


0.5 


5.27 


20 


Si2 


3.4 


22 


0.89 


16 


14 


0.73 


0.45 


4.13 


12 


Si43 


2.0 


31 


0.92 


32 


13.5 


0.61 


0.8 


4.14 


12 


Si6-14 


1.9 


26 


1.22 


26 


15 


0.66 


0.86 


2.35 
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Table [TJ shows that the extracted To and T% agrees rea- 
sonably with h/2ksT. In the theory, the exponents k, and 
C, and the p® c value should be universal. Experimentally, 
there is a trend of k and C" 1 towards smaller values with 
decreasing sample mobility. The trend is rather likely to 
be due to the presence of the long-range components of 
the random potential for the high mobility Si-MOS sam- 
ples. This favors a percolation-type behavior to dominate 
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over the true scaling one [l6j . We mention that variation 
of exponent k with a sample mobility has recently been 
observed in the plateau-plateau transitions of the quan- 
tum Hall regime 17] . Yet another possible reason that 



may complicate the scaling analysis of experimental data 
is an inhomogeneity of the electron density [l8j]. 

In contrast to the one-parameter scaling, in the two- 
parameter case the data should be described by a uni- 
versal scaling function 1Z(X, Y) , which represents a 2D 
surface in the (X, Y, p) space. Using six parameters de- 
termined experimentally, n c , k, C,T 0j1 and p° c (Table 1), 
we calculate the (p/Pc) data versus X and Y. The data 
(64000 points) for all samples collapse close to a sin- 
gle surface in the (X, Y, pj p° c ) space. To model the sur- 
face analytically we use TZ(X,Y) = p° c exp[/i(A)]/ 2 (T) 
and fit the functions fx,2 with five adjustable parameters 
as fx = -X + 0.07^2 + 0.0OIX 3 and / 2 = (1 - Y + 
1.48F 2 )/(1 + 1.88T 2 + 1.72y 3 ). This procedure yields an 
agreement of the model surface with the data within 4% 
rms in the range \X\ < 5 and Y < 3. The X 2 term in 
fx describes the violation of the reflection symmetry for 
large X. The analysis of our p(T) / p c (T) data confirms 
that the reflection symmetry X <-» —X holds within 0.8% 
accuracy in the region \X\ < 0.5, Y < 0.7 (see Fig. 3). 
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FIG. 4: The empirical two-parameter scaling function 
1Z(X, Y) for five samples (dots) and the model surface T = 
XY (grid). The straight line pj p° c =Y at X — represents 
the separatrix. 

The resulting model surface T — XY is plotted in 
Fig. Has a function of X = exp[/ x (X)] and Y = f 2 (Y) 
together with the scaled data p(T,n) / p° c . The data are 
very close to the model surface; the maximal deviation of 
the data from the model surface in the p direction does 
not exceed 0.05p°. Clearly, the surface is curved in the 
both directions. The separatrix p c (T) is represented by 
the straight line on the model surface (Fig. 2]). 

In summary, we have performed a detailed two- 
parameter scaling analysis of the p(T, n) data for five 
Si-MOS samples in the vicinity of the 2D metal-insulator 
transition. The two-parameter scaling function, which we 



have determined experimentally, naturally incorporates 
the temperature dependence of the separatrix p c (T), the 
generalized mirror reflection symmetry with respect to 
p c {T), and the existence of the maximum and inflection 
points for the p(T) dependence on the metallic side. Our 
analysis strongly supports the interpretation of the criti- 
cal behavior observed in the transport data as the mani- 
festation of the quantum phase transition driven by both 
disorder and interaction. The sample to sample varia- 
tions observed in the exponents n and £, and in the crit- 
ical resistivity p° value require a quest for new samples 
in which a random potential will have much shorter cor- 
relation length and the inhomogeneities of the electron 
density will be reduced. 
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